Third-order dissipative hydrodynamics from the entropy principle 
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We review the entropy based derivation of third-order hydrodynamic equations and compare their 
solutions in one-dimensional boost-invariant geometry with calculations by the partonic cascade 
BAMPS. We demonstrate that Grad's approximation, which underlies the derivation of both Israel- 
Stewart and third-order equations, describes the transverse spectra from BAMPS with high accuracy. 
At the same time solutions of third-order equations are much closer to BAMPS results than solutions 
of Israel-Stewart equations. Introducing a resummation scheme for all higher-oder corrections to 
one-dimensional hydrodynamic equation we demonstrate the importance of higher-order terms if 
the Knudsen number is large. 



A causal theory of relativistic dissipative hydrodynamics was first formulated by Israel and Stewart [l[ and has 
been successfully applied to study a wide range of ultra- relativistic heavy- ion collision phenomena @d3j- However, 
recently presented detailed comparisons of solutions of Israel-Stewart equations and kinetic transport calculations have 
^ ■ demonstrated that deviations between them increase with increasing strength of dissipation. Israel-Stewart equations 
can be derived from the divergence of the off-equilibrium entropy current, which is expanded up to second-order in 
dissipative fluxes. It is thus of interest to investigate whether a better agreement between hydrodynamic and kinetic 
transport calculations can be achieved if the entropy current is expanded one order higher than in Israel-Stewart 
theory. 

\Q ' The underlying equation of the kinetic transport theory is the Boltzmann equation 
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which describes the space-time evolution of the phase-space particle distribution f(x,p) — d 3 d J^3 x due to drift and 
diffusion on the left hand side and the collision processes on the right hand side of the equation. A connection between 
O i' the microscopic kinetic transport theory and the macroscopic theory of hydrodynamics can be established using the 
1} , Grad's method, in which the off-equilibrium distribution is approximated by 

foff-eq = fo (1 + e + CM p" + e^pV) (2) 

where fo denotes the isotropic (equilibrium) distribution and the fields e, e M and e M „ are related to the dissipative 
fluxes II, g M and ir^ v . The exact form of e, e M and e M „ can be obtained from the definitions of dissipative currents and 
' the matching conditions 



u^u,AT^ - T£) = 0, u v {N v - N") = (3) 



■ where Tjf" and N" q are the energy-momentum tensor and the particle current in equilibrium. 

In the following we consider a massless gas of Boltzmann particles (gluons) undergoing a boost-invariant one 
dimensional expansion [f|. For the considered system the equation of state is e = 3p . Bulk pressure and heat flux 



vanish identically and the local rest frame off-equilibrium distribution becomes 



e = 6p . 



f(x,p) = d g e~ E ' T ( 1 + — W ^ ) (4) 



with the energy density e and the degeneracy factor for gluons, d g = 16. The matching conditions |(3) | allow to define 
the temperature T for an off-equilibrated system by matching its energy and particle densities e and n to a fictitious 
equilibrium state: 

e = c e q, n = n e q (5) 
For the Boltzmann gas considered here the temperature T in Eq |(4)| is then 

r = f. (6) 
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FIG. 1: Ratio of transverse particle distribution from BAMPS with r\j s = 0.4 to the one calculated by Eq (4) at different times 
with e, T, Tr**" extracted from BAMPS. 



Grad's approximation, Eqs |(2)| resp. |(4)| is essential for derivations of hydrodynamic equations from the Boltzmann 
equation [l5L [20j| resp. from the entropy principle @, It is thus important to quantify how accurate the approxima- 
tion in Eq ](4)| can reproduce the off-equilibrium distribution obtained from the numerical solution of the Boltzmann 
Equation in the partonic cascade BAMPS (T^J. BAMPS has recently been applied to investigate a wide range of 
phenomena such like the buildup of the elliptic flow [l8| , the energy loss of high energy gluon s I19l| , the extraction of 
the second-order viscosity coefficient [9j], and the formation and propagation of shock waves [14| in ultra- relativistic 
heavy-ion collisions. 

For the results presented in this work we use BAMPS calculations with elastic isotropic cross section adjusted in 
such a way that the shear viscosity to entropy density ratio is constant throughout the evolution, as introduced in 
[1, El- The initial condition is a Boltzmann distribution with To = 500MeV at the initial time tq = 0.4/m/c. All 
results are extracted from the central rapidity bin r\ e [—0.1 : 0.1]. 

To quantify the deviations of the expression in Eq |(4)| from the actual distribution in BAMPS we take the ratio 
^Grad PS °f the transverse distributions calculated as follows: 
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In the latter expression 7r denotes the shear pressure, which is the dissipative correction to the longitudinal pressure: 
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T — T eq are 



For the analytic calculation using Eq |(4)| the values of e and T as well as the components of 

extracted from BAMPS. The deviation Sq^ ps is shown in Fig[T]at different times as function of pr- The deviations 
of Grad's approximation from actual distribution in BAMPS do not exceed 10% for the chosen value of r)/s below 
Pt ~ 3GeV. The good agreement of BAMPS distribution with Grad's approximation observed in Fig[T]might indicate 
that the analytic expression in Eq |(4)| can be applied for derivations of hydrodynamic equations. 

Derivation of the evolution equation for the shear stress tensor has been reported by us recently in 10]. The starting 
point for the derivation is the entropy current which in kinetic theory can be written as 
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Writing the off-equilibrium distribution as / = /o (1 
third order in <f> and obtain by a direct calculation 



(comp. Eqs. |(2)| and |(4)[) we expand the logarithm up to 
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with so = 4n — n In A and fa = j-. Note that truncating the expansion at order 4> 2 we obtain the entropy current 
from the Israel-Stewart theory [l], ((| . 

The original Israel-Stewart equations explicitly satisfy the second law of thermodynamics since they are obtained 
directly from the requirement of non-negativeness of the divergence of the entropy current, <9 M s M > 0. Using the same 
argumentation with a third-order entropy current in Eq. |(10)| i.e. calculating its divergence and imposing a linear 
relation between the dissipative flux and the corresponding thermodynamic force, we obtain an extended version of 
Israel-Stewart's equation for tt^ v in which we keep only terms up to third-order in Knudsen number Kn ~ r^d^u^ 
resp. dissipative fluxes [1 01 ] : 

T-n fa n fa ^ \ 2T J 9fa^\T U ) n ° n (n) 

The latter equation constitutes a novel third-order evolution equation for the shear tensor. The notation tt denotes 
derivative with respect to the proper time r. r w is the relaxation time and is the same as in Israel- Stewart theory: 

r T = 2r)fa. (12) 



Neglecting the last two terms in Eq. (11) the second-order Israel-Stewart equation is recovered. To second order the 
equation we obtain does not contain all terms found in recent works 0, EH, EH EH| ■ It is an interesting task for the 
future to understand the differences between various formulations of dissipative hydrodynamic equations, as discussed 
for example in [15j . We have to stress that in our derivations, as presented here and in more detail in [10j, we do not 
use the Boltzmann Equation directly. The terms we obtain are the full set of terms which can be obtained from the 



entropy principle if Grad's approximation, given by Eq. (4) is used. 



For a one-dimensional system with boost-invariance Eq. (4) takes the following form 



* = -^ + A£-3^. (13) 
7v 3 r 27 r er 

In the latter equation 7r denotes the shear pressure which reduces the longitudinal pressure pl- 

Pl = T 33 = p - tt. (14) 

Again, the Israel-Stewart equation is recovered from Eq. |(13)| if the last term is neg lected. Eq. |(13)| has to be 
solved together with the evolution equations for the energy and particle densities. The former is obtained from the 
conservation of the energy-momentum tensor component, d w T v0 = 0. The latter we obtain assuming conservation of 
the particle number current, <9 M -/V M = 0, i.e. assuming a medium in which net particle number is constant. In the 
geometry chosen here, the evolution equation for the energy and particle densities read 

4 e tt ii , . 

e =--- + - , n= — . (15) 

ATT T 

Before discussing the solutions of Eq. |(13)| resp. of the Israel-Stewart's equations we would like to discuss the effect 
of higher than 3rd order contributions to Eq. |(13)[ This analysis will be presented here for a one-dimensional system. 



In order to include all orders of corrections into Eq. (13) we assume they all have the form x n (— ) ^ with n > 3. 
Note that the second and third-order terms in Eq. |(13) are already of this form. Thus, an ansatz for an equation 



containing all orders of corrections can be written in the following way 

7r 4 7r 8 e ^— v /ir\ n e 
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The coefficient x is supposed to be an unknown function of time. Since the equation we consider is supposed to include 
all orders of corrections, it should be applicable in the free-streaming, i.e. r^ — > oo, limit. In the free-streaming limit 
the solutions for e, tt and n are known. Since the gas is streaming free, longitudinal pressure cannot be built up, i.e. 
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FIG. 2: Time evolution of pressure isotropy. The shear pressure is calculated by Eqs. (15) and (16) using \ — (Israel-Stewart) 



X = —3 (extension of Israel-Stewart's equation to third order) and % — ~~ § (approximate inclusion of all order of corrections). 



Pl = 0, which means w = p = |. Using this, one finds that the energy density evolves according to Eq. |(15)| as 
e(r) = cqTq/t. Using these solutions in the tv — > oo limit in Eq. (16) one obtains 



(17) 



This value is the result of resummation of higher-order terms in the heuristic ansatz Eq. (16) The equation obtained 



this way includes all orders of corrections, but only approximately. If only third-order terms are included, one obtains 
X = 3, which corresponds to Eq. (13)| The Israel- Stewart equation is obtained by setting \ = 0. 



In the following we present the solutions of hydrodynamic equation of second (Eq. 



orders (Eq. (1 



(16) with x = 0) and third 
with x — —3) and of the approximation of all orders (Eq. (16) with x = § )■ The observable we use 



to quantify the deviations from equilibrium is the pressure isotropy 



The results are compared to BAMPS 
For the comparisons presented 
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calculations using thermal initial conditions with Tq — O.bGeV and To = 0.4/m/c. 
here only elastic processes with isotropic cross section are included in BAMPS. The cross section 022 is parametrized 
to keep the rj/s value constant, as has been already done in HoL [Til] : 
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where 4n — n In A = s. The method of derivation of hydrodynamic equations which has been presented here does not 
allow to obtain an analytic expression for the shear viscosity coefficient. We thus rely on the expression for 77 obtained 
in [2l| which is strictly speaking valid for the Israel- Stewart, i.e. second-order, theory. 

Time evolution of the pressure isotropy is presented in Figj2j It was demonstrated in [l4[ that the evolution of 

the system is governed by the Knudsen number Kn = t^/t rs 6^^. Since T w T (t /t) 3 , the Knudsen number 
depends in the situation considered here on tq,t, T and rj/s. Since all results are obtained using the same initial 
conditions, the evolution depends only on the chosen value of rj/s, which is thus the direct measure of the strength 
of dissipative effects. The deviations between the Israel-Stewart and BAMPS results become considerable already at 
rj/s = 0.4. At rj/s — 3 the Israel-Stewart's equations (x = in Eq. |(16)[ ) lead to a negative longitudinal pressure 
which is not physical for the considered setup (This phenomenon has been as well investigated in Q)- The absolute 
deviations are reduced if one extends Israel-Stewart's equations to third-order, Eq. |(13) resp. x = — 3 in Eq. |(16)| 
Moreover, the negative longitudinal pressure does not occur in the solutions of the third-order equations. Especially 
at late times, where the relaxation towards equilibrium sets in, the third-order results on pressure isotropy are in very 
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good agreement with BAMPS solutions. At early times the third-oder and BAMPS results still deviate since there 
the expansion scalar d^u^ = 1/r, is still large and thus all orders of corrections have to be taken into account. Indeed, 



at early times solution of the approximate all-order equation, \ — — § in Eq. (16) is in very good agreement with 



BAMPS results. This is due to the fact that the value x — ~ § nas been obtained from the — > oo or alternatively 
Kn — > oo limit, to which the system is close at early times. 

In this study we have presented an extension of the Israel and Stewart's entropy based approach to third order in 
dissipative fluxes. For this study we have considered a one-dimensional boost-invariant gas of massless Boltzmann 
particles (gluons) . Results of hydrodynamic calculations have been compared to solutions of the Boltzmann Equation 
from the partonic cascade BAMPS. Although the Grad's approximation, upon which the derivation of Israel-Stewart's 
equation is based, has been shown to describe transverse spectra from BAMPS with remarkable accuracy, the solu- 
tions of Israel-Stewart equations demonstrate large deviations from BAMPS. Inclusion of third-order terms into the 
evolution equation for shear tensor reduces the deviations between hydrodynamic and BAMPS results considerably. 
Up to r]/s = 0.4 the third-order solutions are in very good agreement with BAMPS. In order to estimate the effect 
of all orders of correction, we have introduced a resummation scheme which allows to represent the infinite series of 
higher-order terms by one single term. The solution of this equation is in very good agreement with BAMPS results 
even at rj/ s = 3, when the Knudsen number is very large, which underlines the importance of higher-order corrections 
at early times of evolution. It is still an important task for future studies to understand the differences between the 
equations obtained here from the entropy principle and the second-order equations presented in recent publications 
by different authors [!, [ID, [TH, [2(| . It is as well important to investigate whether a possible extension of Grad's 
approximation can further improve the agreement between hydrodynamic and kinetic transport calculations. 
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